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I. Introduction 



The steady aerodynamic characteristics of the jet flap and its deriva- 
tives have been studied rather extensively over the past two decades and a 
considerable body of information including a general computer program (Ref, 12) 
for its prediction is now available. The unsteady jet flap, characteristics are 
much less understood although the potential of the jet flap as a fast acting 
lift control device was recognized quite early. In the last few years there 
was a noticeable increase in activity on this problem ranging from experimen- 
tal/theoretical studies of the jet flap's frequency response to applications 
for wind tunnel dynamic response investigations as well as aircraft and 
helicopter gust load and vibration control. It is the objective of this lec- 
ture to present an up-to-date review of these activities. 
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2. Unsteady Jet Flap Theory 

2.1 Basic Theory 

Assume an airfoil that has a jet emerging at the trailing-edge . Both the 
airfoil and jet are free to execute a small amplitude but otherwise arbitrary 
time-dependent motion. The jet is assumed to extend infinitely far downstream 
and both the main-stream flow and the flow in the jet are taken to be inviscid, 
incompressible and irrotational. Furthermore, it is assumed that the local jet 
velocity v is very much greater than the local velocity of the main stream u 
in the vicinity of the jet. The total pressure in the jet is, for v>>u , 
greater than that in the main stream which requires in the absence of viscous 
mixing and because of the continuity of static pressure across the jet bounda- 
ries that a vortex sheet satisfies the velocity discontinuity there. 

In the main stream the unsteady Bernoulli equation must hold 
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where p, (j>, u are the local static pressure, velocity potential and velocity, 
respectively, is the constant density, the subscript zero denotes conditions 

at infinity. Similarly, in the jet one has 
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where p, (J)^, v are the local static pressure, velocity potential and velocity 

in the jet, p, is the constant but in general different jet density, p and 
*J o 

are the static pressure and velocity at infinity downstream independent of 
time. This is based on the assumption that sufficiently far downstream at a 
given instant of time after initiation of the motion, the jet must return to 
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its initial undisturbed position, requiring also that its slope and curvature 
vanish. Therefore, one has the condition 
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The following argument is advanced to justify this statement (Erickson 1962): 
Since the jet is Impermeable . and extends to infinity downstream a particle 
above (or below) the jet initially must always remain above (or below) . If 
the jet were displaced from its initial position at infinity downstream, an 
infinite amount of work would have been done in a finite time to move the 
infinite amount of air above and below the jet. 

Also, since the flow in the jet is irrotational one has 




where r is the radius of curvature of the streamlines. Integration of this 
equation then gives 



rv = Constant (4) 

To obtain the pressure difference across the jet consider an incremental 
jet element shown in Fig. 1. Subscripts 1 and 2 refer to upper and lower jet 
boundaries. R(s,t) is the radius of curvature, 6(s,t) the jet thickness at 
position s and time t. 
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From Eq* (2) one obtains 
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However, since the radii of curvature are normal to the streamlines inside 
the jet they are equipotentials of the jet flow. Hence 



<|) (s,t) - (f> (s,t) = 0 

Jl J2 



Defining the mean jet velocity as 



V(s,t) = 



Vji^Csjt) 



the irrotationality condition Eq. (4) becomes 
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and therefore 
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Substituting the above equations back into the equation for the pressure 
difference across the jet one finally obtains 

PjV^(s,t) 5(s,t) 

Pj(s.t) - PjCs.t) 
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where the jet momentum is 



J(s,t) = PjV (s,t) 6(s,t) 



( 6 ) 



Eq. (5) shows that the pressure difference across the jet is proportional to 
the jet momentum flux and inversely proportional to the radius of curvature of 
the jet. This is analogous to the steady jet flap result and was achieved by 
the assumption v»u causing the vanishing of the velocity potential difference 
across the jet. 

Having expressed the pressure difference across the jet it is next re- 
quired to determine the jet vortex strength. Using Eq. (1) for the main stream 
across the jet one obtains 
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Defining the mean velocity of the main stream across the jet as 

u^ (s,t) + u,(s,t) 

U(s,t) ^ 

the velocity difference can be expressed from the above equation together with 
Eq. (5) as 
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In contrast to the steady case one now has an additional term in the form of 
the time derivative of the velocity potential difference. 

For the limiting case of an infinitely thin, high-speed jet whose jet 
thickness 5 -v 0 tends to zero it is possible to relate the jet vortex 
strength to the velocity difference across the jet, hence 

Yj(s,t) = u^(s,t) - U 2 (s,t) (8) 



For a rigorous proof compare Erickson (1962). In this limit V(s,t) must 
become infinite for finite momentum flux although the flow is still considered 
incompressible. Eqs. (7) and (8) then become 
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To complete the problem formulation it is necessary to enforce the flow 
tangency (kinematic) boundary condition on the airfoil and jet. The linearized 
approximation of this condition states that the downwash is given by the 
linearized convective derivative of the airfoil and jet ordinate (see e.g. 
Kararacheti 1966), i.e. 
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Where the airfoil extends from 0<x<“. 
Linearizing also the jet curvature 
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one obtains for Eqs. (5) and (9) 
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and 



Yj(x,t) = 



_ 1 IT p ^ 3 Y(x,t) _ JL 9A<t>(x,t ) 

2 o J .2 U 9t 

OX O 



C<X<oo 



(13) 



Where s has been replaced by x due to the linerization and the jet momen- 
tum coefficient 
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has been introduced. 
Using 
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the potential difference can be expressed from Eq. (13) as 
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and differentiating this equation with respect to x gives 
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Representing now both the airfoil and jet by vortex distributions the 
downwash due to these vortex distributions is found from the Biot-Savart law 
giving the well known equation 
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w(x,t) = - 



1_ 

2 tt 



00 

^ Y(g,t) dg 

C - X 



0<X<oo 



(17) 



If one writes 
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then Eq. (17) can be inverted using one of several techniques (see Heaslet and 
Lomax 1954, Carleman 1922 or Cheng and Rott 1954) leading to the following 
equation on the semi-infinite interval x>c 
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w (x,t) is known from the flow tangency on the airfoil Eq. (10) for known 
w 

airfoil motion. Along the jet this condition is 
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where Y(x,t) = hj(x,t) is the jet deflection below the x-axis. 
one has the dynamic jet condition Eq. (16) 



In addition. 
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Eqs. (18), (19) and (20) are the basic starting equations for the solution of 
the unsteady jet flap problem within the linearized thin-jet approximation. 
They are three equations for the three unknowns Yj(x,t) , Wj(x,t) and hj(x,t) , 

but with prescribed wing motion w (x,t) and prescribed variation of 

w 

the jet slope at the airfoil trailing-edge. For a more detailed derivation 
and discussion of this problem formulation compare Spence 1961 and Erickson 
1962. 

2.2 Spence's Oscillatory Flow Solution 

Assuming purely harmonic time dependence Spence (1965) has given an 
approximate solution of the above system of equations. Three cases are of 
interest, namely a) the stationary airfoil with oscillating jet flap, b) the 
periodically plunging airfoil c) the periodically pitching airfoilo 

For case a) the airfoil remains fixed at zero angle of attack, hence 






h (x,t) = 0 = w (x,t) 
w w 



( 22 ) 



In cases b) and c) the airfoil moves but the jet emerges at all times 
tangentially at the trailing-edge. 
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Replacing hj(x,t) and 6j(x,t) by e"*" ^o 



Eqs, (18) and (19) combine into 
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Eq. (20 becomes 
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and the boundary condition is 



hj(c) = T = Const, hj(c) = 0 

Spence proposed a ”weak-jet" approximation to this pair of equations by assuming 
the jet strength parameter y to be small with respect to unity or, equiva- 
lently, Cj to be small with respect to 4. "Stretching" the trailing-edge 
region by the coordinate transformation 

X = (x-c)/yc (25) 

he derived an approximation to Eqs. (23) and (24) which is dependent only on 
the single parameter 

V = yw (26) 

The resulting singular third-order intcgro-differential equation then becomes 
mathematically tractable. Spence found solutions to exist only for values of 
V ^ 2 and computed jet shapes and lift forces for both the oscillatory jet case 
a) and the plunging case b) for a range of values of v, Erickson (1970) later 
extended this approach to the harmonically pitching airfoil, case c) , and 
presented lift and pitching moment results. 

According to Erickson it is important to note the "boundary layer" 
character of this solution which was obtained by the stretching transformation 
Eq. (25) and then setting y equal to zero whenever it appeared explicitly in 
the transformed equations. The neglect of the finite y terms, together with 
the transformation, emphasizes the jet just downstream of the trailing-edge 
and restricts the validity of the solution to high frequencies. This follows 
from the assumption that the characteristic length c and the characteristic 
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wavelength of the oscillation, given by 2irc/ , must both be of the same 

(jj 

order. This implies that w be of order 2ir or greater. 

An interesting finding is the breakdown of the solution for values of 
V > 2, Although no completely satisfactory explanation seems obvious, Spence 
hypothesizes the existence of a critical "Strouhal number" indicating a wake 
instability, 

2,3 Potter's Approach 

Potter (1972) approached the problem as an initial value problem in time 
using point vortex distributions and thus reducing Eq, (17) to finite difference 
form. The jet is initially specified as a straight horizontal line with zero 
vorticity, after which time the jet exit angle is either instantaneously 
deflected or continuously oscillated. By accounting for all the changes in 
bound airfoil vorticity and by assuming that the shed vortices are propagated 
domstream with the free-stream velocity, conservation of circulation is main- 
tained for each time step and the jet shape is determined by computing the 
downwash and the resultant velocities at each jet vortex point. As is common 
with point vortex methods, erratic vortex motions were found to build up 
with increasing time steps. Although no completely satisfactory solution to 
this problem was found, the adoption of an averaging scheme over several 
neighboring vortices gave seemingly satisfactory results for jet momentum 
coefficients smaller than 0,1, Above this value the procedure was unstable. 
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3. Unsteady Jet Flap Experiments 

In the late 1960 's several investigators started to conduct wind tunnel 
tests on oscillating jet-flapped airfoils in an effort to determine their 
unsteady aerodynamic characteristics. 

3.1 The Work of Simmons and Platzer 

Simmons and Platzer (1970) studied a two-dimensional NACA 0012 airfoil 
of 12-inch chord and 18-inch spano The aft 15% of the chord was removed and 
replaced across the full span by a circular steel tube with 0.5-in. o.d. as 
shown in Fig. 2. There was no gap between the tube and the airfoil and the 
tube was free to rotate in bearings at each end. It served as a plenum for 
the jet flap which was formed by a row of 0.029-in. diameter holes spaced 
0.25 inches apart along the trailing edge of the tube across the full span. 

Fig. 3 shovJs the test set-up. The airfoil was clamped between two rectangular 
end plates extending from the tunnel roof to the floor and extending 0.6 
chord lengths upstream and 2.6 chord lengths downstream of the airfoil. The 
tests were performed in the Lockheed-Georgia Low Speed Wind Tunnel which has a 
test section 30 in. high, 43 in. wide and 58 in. long. The jet flap’s oscilla- 
tion was accomplished by extending the plenvim tube through a bearing in one of 
the end plates and coupling it to the shaft of an electrohydraulic servo-motor 
fastened to the rear wall of the tunnel. The jet deflection angle could then 
be varied, at frequencies up to 30 cps, by rotational oscillation of the plenum 
tube. The airfoil instrumentation consisted of 21 pressure tappings that 
were located on both the upper and lower surface along the chord at mid-span. 
These tappings were connected through a scanning valve to a single pressure 
transducer located outside the tunnel. Thus, the expense and calibration 
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of a large number of transducers located directly on the airfoil surface was 
avoided. Instead, the frequency dependent transfer function of the tube- 
scanning valve system had to be determined to reduce transducer outputs to 
actual pressures on the airfoil surface. This technique was first pioneered 
by H. Bergh (1964) of the National Aerospace Laboratory. 

The tests were performed with the model set at zero angle of attack and 
with the jet flap oscillating about the airfoil chord line at 13 degree ampli- 
tude. The tunnel speeds were 103 and 51 fps giving jet momentum coefficients 
of 0.14 and 0,58. Typical differential pressure distributions are shown in 
Figs, 4, 5 and 6 for three frequencies of oscillations, 1, 10 and 22 cps. These 
pressure distributions were then integrated to give the coefficients and their 
phase relative to the plenum tube position which are presented in Fig. 7. 

3.2 The Work of Trenka 

Trenka (1970) investigated a two-dimensional NACA 0015 airfoil having a 
12-inch span and a 6-inch chord with a 20% flap over the entire span. The 
jet exhausted over the upper flap surface from a span-wise slot at 80% chord. 

The airfoil was bolted to a main support tube which was supported in bearings 
located at 38% chord. By means of electromagnetic shakers and suitable 
shaker arms the wing and flap could be excited so that jet flap oscillations 
as well as wing pitching oscillations could be generated. The wing and 
support tube were attached to a three- component strain-gage balance permitting 
the measurement of the loads perpendicular and parallel to the wing chord 
plane and the pitching moments about the wing pitch axis. The tests were con- 
ducted in the subsonic leg of the Cornell Aero Lab High-Speed Wind Tunnel 
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having a test section 17 inches wide, 32 inches high and 32 inches long. The 
dynamic tests were performed with and without a tunnel velocity in order to 
provide information about the balance inertia loadings at zero wind speed. 

Both, wing and flap were set at zero angle and the tests were run with and 
without jet blowing. Also, a number of pressure transducers were mounted on 
the upper and lower airfoil surface to give static and oscillatory pressure 
information. Results are shown in Figs. 8 and 9. 

3.3 The Work of Takeuchi 

A third series of tests was completed at about the same time by Takeuchi 
(1970). Again, a dynamic balance was used to determine the oscillatory forces 
and moments rather than direct pressure measurements. The semi-span wing 

model had a rectangular planform with a 7.5-inch chord and a 24-inch semi- 

span. The profile approximated a NACA 632A015 airfoil. A plain trailing 
edge flap of 0.65-inch chord was attached to the wing over the full span for 
the purpose of varying the jet angle which emanated from blowing slots at the 
top and bottom of the wing trailing edge. An eccentric arm of adjustable length 
driven by an electric motor was used to generate flap oscillations of variable 
amplitude. The wing was mounted on a strain gage balance which provided a 
measure of the lift at any instant of time. The tests were performed in the 
3 ft X 3 ft subsonic tunnel of the Pennsylvania State University. The tunnel 
speed was kept constant at 56.5 ft/sec, permitting jet coefficients up to 

0.5. Using galse walls various aspect ratios could be simulated, ranging 

from 2.1, 4.1, 6.4 to infinity. Some of the results are shovm in Figs. 10 and 

11 . 
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3.4 The Work of Kretz 

Kretz (1973) completed another series of tests closely paralleling the 
work of Simmons and Platzer (1970) , The NACA 0012 airfoil tested had a span 
of 1000 ram and a chord of 400 mm. The last 20% were replaced by a flexible 
flap to achieve various jet deflection angles. The measuring technique was 
based on the previously described NLR-technique (Bergh 1964) requiring the 
tube transfer functions. The tests were conducted in the ONERA S2L tunnel 
(no test section dimensions were given) at a wind speed of 40 m/sec, per- 
mitting tests at jet momentum coefficients of 0.25, 0.5 and 1.0. The results 
are shown in Fig. .12. 

3.5 The Work of Simmons 

Very recently, Simmons (1974) completed another series of tests on a 
NACA 0012 airfoil which had a chord of 150.9 mm and a span of 461 mm. This 
experiment was performed in the Low Speed Wind Tunnel of the University of 
New South Wales, Australia, which has essentially a 0.461 mm x 0.461 ram 
square test section except for corner fillets. Again, the measuring technique 
was the same as the one previously used by Simmons and Platzer (1970). The 
airfoil could be excited into sinusoidal pitch or plunge oscillation by a 
precision Scotch yoke mechanism acting through proper linkages, while the 
jet flap was held stationary with respect to the airfoil. The jet momentum 
coefficient could be varied up to 1.76 running at tunnel speeds of 12.3 and 
24.6 m/sec corresponding to Reynolds numbers of 1.29 x 10^ and 2.58 x 10^. 

Some of the results obtained in these tests are shown in Figs. 13 and 14. 
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4. Comparison and Evaluation 
of Tests and Theories 

4.1 Fixed Airfoil with Oscillating Jet Flap 

The results of Spence, Simmons and Platzer, Takeuchi, Potter and Kretz 
are compared in Figs. 15 and 16. It is immediately apparent that there is 
relatively good agreement between Spence's theory and Takeuchi’s experiment 
on the one hand, and between Potter's theory and the experiments of Simmons 
and Platzer and Kretz on the other hand. The large discrepancy between these 
two sets of theories and experiments must remain largely unexplained until 
further experimental and theoretical work is accomplished. At this time only 
the following observations can be offered: 

a) As previously explained Spence's theory appears to be valid only for 
reduced frequencies co > 2ir and a comparison of this theory in the 
lower frequency range may therefore be misleading. 

b) The deviations between the experiments of Simmons and Platzer and 
Kretz may be caused by the different trailing-edge and jet configura- 
tions. Simmons and Platzer used a blunt trailing-edge and many 
discrete holes rather than a continuous slot for jet blowing, whereas 
Kretz had a sharp trailing-edge and a continuous jet sheet. 

c) The major unresolved discrepancy remains between the dynamic balance 
measurements of Takeuchi and the pressure measurements of Simmons and 
Platzer and Kretz which produced quite opposite trends for both lift 
amplitude and phase angle with frequency. 
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4.2 Periodically Pitching and Plunging Jet Flapped Airfoil 



Simmons (1975) compared his recent results (section 3.5) with Spence's 
predictions and Trenka's experiments. Again, the theory disagrees with the 
experimental trends over the reduced frequency range («<27 t) tested but this may be 
explainable by the inapplicability of the theory in this range. However, the 
irreducibility of Spence's results for lift and pitching moment to the well- 
known jet-off results as the jet momentum coefficient approaches zero is 
particularly disturbing and requires further investigation. Trenka's experi- 
mental results, on the other hand, are too limited to allow meaningful 
comparisons. 
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5. Applications 



The potential of the jet flap as a fast acting lift control device seems 
to have been first recognized by W. R. Sears in the late fifties (according to 
a remark by Spence in Ref. 19) . Since then a number of applications have been 
made or proposed which we will briefly discuss. 

5.1 Air Stream Oscillator 

The need for practical and efficient airstream oscillators in wind tunnel 
dynamic response studies is well recognized. Several devices have been investi- 
gated in recent years to find alternatives to the presently used oscillating 
vane arrangements (such as the one used in the Langley Transonic Dynamics 
Tunnel, Ref. 1 ) Jacobs and Platzer (1969) proposed mechanically or fluidically 
oscillated jet flaps for airstream oscillation. A typical configuration is 
shown in Fig. 17 which was tested by Simmons and Platzer (1970). Further work 
along these lines was recently pursued by Viets (1975) and Ham (1975) . Viets 
(1974) developed a new fluidic oscillator V7hich appears to be well suited for 
high frequency gust generation. Ham (1975) adopted the circulation control 
principle rather than the jet flap principle. His set-up is similar to the one 
used by Simmons and Platzer but has the advantage of significantly lower re- 
quired blowing rates. Test data on the performance of these devices can be 
found in Refs. 17 and 18. 

5.2 Gust Response and Vibration Controller 

In recent years extensive studies have been conducted on active load 
alleviation, mode stabilization and vibration control of both aircraft wings 
and helicopter rotor blades using conventional control surfaces and spoilers 
as control devices (see e.g. Refs. 13,20). Helicopter blades can experience a 
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particularly complex non-uniforra flow field due to the interception of vortices 
produced by preceding blades. The forces produced by these flow non-uniformities 
may excite higher modes of harmonic vibration and high frequency stresses 
leading to premature blade failure. Closed-looped control techniques utilizing 
pressure, lift or strain sensing on the blade section therefore appear as an 
attractive method of vibration and load control. The application of the jet flap 
to this problem has been studied extensively by Giravions Dorand (Refs. 11, 13) 
and to a lesser extent by Honeywell Inc. (Ref. 2). Giravions Dorand con- 

ducted a series of wind tunnel tests on a jet flap rotor to reduce the vibrations 
and stresses by multi-cyclic jet flap control. At a tip speed ratio of 0.4 
reductions up to 50% were attained. Concurrent theoretical studies showed that 
further reductions up to 80% are possible. Honey\\^ell investigated the Variable 
Deflection Thruster (VDT) concept (Fig, 18) . The VDT allows control of both 
the angle and magnitude of a jet by mixing of two separate jets. The two jets 
are supplied from two separate plenum chambers and meet on the surface of the 
blunt trailing-edge at a position determined by their momentum difference. A 
different jet actuation scheme was developed by Simmons and Platzer (1970). 

The effectiveness of this device to alleviate the gust response of an airfoil 
having the pitch degree of freedom only was demonstrated by Platzer, Deal and 
Johnson (1975). 
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Plgure 1 - Flov in an Element of the Jet (Schematic) 
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Fig. 2 Airfoil Used by Simmons and Platzer (Ref. 18) 




Fig. 3 Jet Flap Model Installation in Wind Tunnel (Ref. 18) 
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4 Differential Pressure Distributions at 1 c.p.s, 
Obtained by Simmons and Platzer (Ref. 18) 
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Differential Pressure Distributions at 10 c.p.s, 
Obtained by Simmons and Platzer (Ref. 18) 
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Differential Pressure Distributions at 22 c.p.s. 
Obtained by Simmons and Platzer (Ref. 18) 
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Fig. 7 Measured Lift Amplitude and Phase Angle 
Obtained by Simmons and Platzer (Ref. 18) 
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, PHASE ANGLE, deg 




Fig. 9 plot of pitching MOMENT PHASE ANGLE VS REDUCED FREQUENCY 
From Trenka and Erickson (Ref. 22) 
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From Takeuchi (Ref. 21) 



and «- 



28 




Fig. 11 Variation cf Phase Angle with Reduced Frequency 
at Constant c, of 0.1 and 0.42 AR = 



From Takeuchi (Ref. 21) 
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Fig. 12 Measured Lift Amplitude and Phase 

Obtained by Kretz (Girarions-Dorand, Ref. 11) 
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Fig. 13 In-phase and Quadrate Components of 

Differential Pressure Distribution for 

Oscillatory Pitching 

Obtained by J. M. Simmons (Ref. 16) 



31 



o« 



10 



8 



dCi 

8o^ 



o 



l^o o 

X 

I X X 



4 

0° 



- 10 ' 






- 20 ' 



-30' 



O 



Cj = 1-75, V = 12-3 m s 



-1 



O Cj = 0-435. V = 24-5 m s 
X Cj = 0-215. V = 24.5 m s 



□ Ci = 0 



V = 24 • 6 m s 
o 



0-0 



0-2 0-4 0-6 0-8 1-0 

DiMENSIONLESS FREQUENCY - k 



o 



□ 

□ 

3 


X 

□ 

! 


O 

X 

□ 

1 


o 

X 

□ 

♦ 


1 J 


1 

1 

1 __ _ 




1 


i 


i 


' ' ' o 












o 

— 








o 








O 






O 


X 








! 


J 


□ 

_J 


J ^ ! 



1-2 



Fig. 14 Measured Lift Response for Oscillatory 
Pitching of Airfoil 

Obtained by J. M. Simmons, to be published 
in AIAA-Journal 
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LIFT AMPLITUDE RATIO 



UNSTEADY LIFT MAGNITUDE 
OSCILLATING 2-D JET FLAP 




Fig. 15 Comparison of Results for Lift Response 
of Airfoil with Oscillating Jet Flap 
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UNSTEADY LIFT PHASE ANGLE 
OSCILLATING 2-D JET FLAP 




REDUCED FREQUENCE, K 

Fig. 16 Comparison of Results for Lift Response 
of Airfoil with Oscillating Jet Flap 
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Fig. 17 Gust Generator Used by Simmons and Platzer (Ref. 18) 
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Air-foil Section with a Variable Deflection 
Thruster Trailing Edge 




Varialrk' Deflection 'rhrustei- Schematic Diagram 



Fig. 18 From Bailey and Hammer (Ref. 2) 
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